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ABSTRACT 


This paper describes a unified theory of design sensitivity analysis of linear 
and nonlinear structures for shape, nonshape and material selection problems. The 
concepts of reference volume and adjoint structure are used to develop the unified 
viewpoint. A general formula for design sensitivity analysis is derived. Simple 
analytical linear and nonlinear examples are used to interpret various terms of the 
formula and demonstrate its use. 


1 . INTRODUCTION 


Design sensitivity analysis gives trend information that can be used in the 
conventional or optimal design process. The subject, therefore, has received 

considerable attention in recent years. For a thorough review of the subject Refs. 
1 and 2 should be consulted. 

The present paper describes a unified variational theory of design sensitivity 
analysis of linear and nonlinear structures (geometric as well as physical non- 
linearities) including shape, nonshape and material selection problems. The adjoint 
variable approach is utilized although the direct differentiation method can be also 
easily developed. In Section 2, equations of continuum mechanics for nonlinear 
analysis are summarized. They are needed in design sensitivity analysis. A unified 
viewpoint for shape and nonshape design sensitivity analysis is described in Section 
3. The concept of a reference volume is explained in Section 4. The variational 
theory of design sensitivity analysis using adjoint variable approach is developed 
in Section 5. The theory is used to solve several simple analytical problems in 
Section 6. Finally concluding remarks are given in Section 7. 


2. NONL INEAR ANALYSIS 


Nonlinearities in structural systems can be due to large displacements, large 
strains, material behavior and boundary conditions. Consistent theories to treat 

these nonlinearities have been developed^’ 1 *. We will use the developments and 
notations of Ref. 4, and follow the Total Lagrangian (or Lagrangian) formulation, 
although updated Lagrangian formulation can also be used. One of the major diffi- 
culties in describing nonlinear analysis is the complexity of notation. We will 
mostly use standard symbols from the literature for various quantities. Matrix and 
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tensor notations will be used. One major departure from linear analysis is that 

quantities must be measured in a deformed configuration. Also, a reference con- 
figuration for the quantities must be defined. We will use a left superscript to 
indicate the configuration in which the quantity occurs and a left subscript to 
indicate the reference configuration. 

A starting point for theory of nonlinear analysis is the principle of virtual 
work for the body in the deformed configuration at time t (load level t): 


Io v ° dv - o r - stu ° dv ■ Jo. 0 T - 5 ‘“ \ ■ 0 


'°v 0 


(1) 


where left subscript 0 refers to the undeformed configuration, a refers to the 

standard tensor product and 

= undeformed volume of the body 

qS = Second Piola-Kirchhof f stress tensor 

qE = Green-Lagrange strain tensor 

= body force per unit volume 

= displacement field 

qT = surface traction specified on part of the surface r^, 

°r = surface of the body 
6 = variation in the state fields 


Let u° be the specified displacement on the part r of the surface. The variations 

of the state fields in Eq. (1) are arbitrary but kinematically admissible. They can 
be replaced by any kinematically admissible fields. In particular they will be 
replaced by adjoint structure state fields in later derivations. The virtual work 
equation can also be written using Cauchy stress tensor and other quantities 
referred to the deformed configuration. Transformation can be used to recover 
Cauchy stresses from second Piola-Kirchhof f stresses and vice versa. However, in all 
the derivations given in this paper we will use the undeformed configuration as the 
reference conf iguration . 


The Green-Lagrange strain tensor is given as 


t 

0 


e 


• Ho'*- 1 * 


r o T 

( 0 V » 



(2) 


The nonlinear stress strain law, in general, can be written as 

qS = ♦(Je,b) ( 3 ) 

where b is a design variable. Note that for many applications, functional form 
for $ is not known. In numerical implementations, the explicit form is not 
needed. Only an incremental stress-strain relation is required. 
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There are 


Equations (1) to (3) are nonlinear in the displacement field fc u. 

several methods for solving such system of equations. ^ The incremental/iterati ve 
procedure based on Newton methods is the most commonly used and effective procedure. 
This will be summarized here. In the derivation of the procedure, it is assumed 
that equilibrium is known at t and it is desired at t+At. The state fields are 

ii 

decomposed as 


t+At t ^ t+At_ t 

u ■ “ ♦ “• 0 s ■ o s o Sl 


o s ■ * >E o £ 


t+At t . ^ t+At, t t+At t 

o c ■ o e o e * o f o f o f; o T 0 T * 0 T 


(«) 


where u = increment in the displacement field 

qS = increment in the Second Piola-Kirchhof f stress 

Q e = increment in the Green-Lagrange strain 

= increment in the body force 

qT = increment in the surface traction 

Variation of the strain field is given as 


.t+At 

s o e ■ V 

The incremental strain field from Eq. (2) is given as 

o e = o e + o n 


0 2 0 


I[ A Vu T + ( n Vu T ) T +( n Vu T ) (yu T ) T +( n vV) LVuY] 


(5) 


(6) 

(7) 


" 2 [( 0 VU 



( 8 ) 


Substituting Eqs. (4) - (6) in the virtual work principle, Eq. (1), written at t+At 
and using the fact that state at t is in equilibrium, we obtain the following incre- 
mental virtual work principle: 

Jfr + o S K e ° dV " fo f - 6U ° dV ' fo T - 6U ° dr T = ° (9) 


Equation (9) is still a nonlinear in incremental displacement field u. It is 
linearized by assuming 


6 0 E " 6 o e; 


0 S = *,e‘0 e 


( 10 ) 


and iteration is used within the load increment to satisfy the equilibrium exactly 
at t+At. The finite element procedure has been used to implement the preceding 

equations into a computer program ADINA.^ 
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3. UNIFICATION OF DIMENSIONAL AND 
SHAPE DESIGN SENSITIVITY ANALYSIS 


In the literature, shape and dimension design sensitivity analysis problems 
have been treated independently. In the shape problem, domain of the problem is 
allowed to vary whereas in the dimensional problem domain is fixed but cross- 
sectional dimensions are allowed to vary. It will be seen here that when varia- 
tional formulation is used and volume integrals are used, there is no distinction 
between the two problems. 

Consider the general functional requiring design sensitivity analysis: 

* ' k(b)» S '° C ' tu ’ b) ° dV * ^° r (M° r ’ b> ° dru * ^ S(t “’ b)0dr T <") 

r u b r T (b ) 

It can be seen that when design b is changed, the volume of the body as well as its 
surface change. As examples, consider optimal design of two simple bodies shown in 
Fig. 1. Are these shape or dimensional optimization problems? Our contention is 
that although length of the members is not treated as a design variable in these 
problems, volume of the body changes whenever any of the indicated design variables 
is changed. We must account for variations of the domain of the body while writing 
variations of the functional iJj in Eq. (11). Thus the variational concept for design 
sensitivity analysis is slightly different from the corresponding concept used in 
purely analysis problems where domain of the body remains fixed (at least in linear 
problems). This distinction is important in maintaining generality of the varia- 
tional design sensitivity analysis theory where variation of the domain should be 
always considered. 



DESIGN VARIABLE : A(x) 



DESIGN VARIABLES: w(x), d(x) 


Figure 1. Examples of Optimal Design 
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4. CONCEPT OF REFERENCE VOL WE 


The concept of a reference volume is extremely useful in problems where the 
volume of the body is changing. The idea, introduced recently in Ref. 7, is to map 
volume of the body in various configurations to a reference volume V. This is shown 

in Fig. 2. The original volume of the body °V(b°) moves to a volume t V(b°) under a 
nonlinear motion. However, both the volumes can be mapped to the fixed reference 

volume V under the mappings F^b 0 ) and F ? (b°) respectively. The design process 

changes shape of the body so that its volume becomes ^V(b^) at the new design b^ . 

f 1 

This volume moves to V(b ) under the nonlinear motion. Both these volumes can also 
be mapped to the fixed reference volume V. 

The concept of reference volume is also quite useful in design sensitivity 
analysis. All the integrals of the problem are transformed to the reference volume 
using the proper transformation of the independent variables. The mapping to the 
fixed volume keeps changing under state or design variations. However, the refer- 
ence volume never changes. Thus, when variations of various integrals are taken, 
the variations of the reference volume need not be considered. In numerical imple- 
mentations, this concept is also very useful. It allows us to discretize the design 
problem into design elements that keep the same shape even when the real shape for 



117 




the structure changes during the optimization process. Using the transformation of 
independent variables, various expressions are given as 


Virtual Work Equation at Load Level t: 

JqS.SqE J dV - JJf .6 fc u J dV - JJt.sSi 3 dr T = 0 
Incremental Virtual Work Equation at Load Level t+At: 
J(Js+ 0 S).6 0 e J dV - J Q f.6uJ dV - J q T.6u J dr T = 0 
Green-Lagrange Strain Tensor: 

q e = ^[x^vV) + ( p vV) T x + x T ( vV) ( r vV) T x] 


Incremental Strains: 
IroT, „T 


Q e - + ( r Vu T ) T X ♦ X T ( r Vu T > ( r vV) T X 


,t T J; 


0 n = |[x T ( r Vu r ) ( r Vu T ) T x] 


+ X T ( p Vu T ) ( r VV) T x] 


T.T-- 


( 12 ) 


( 13 ) 


(U) 


(15) 

( 16 ) 


Functional for Sensitivity Analysis: 

* = / GfJs.Jje.V) J dV + J g(jjT,b) 3 d? u + J h( fc u,b) 3 df T (17) 

Jacobian of Transf ormation: 


~,0 0 0 , 

3( x , y, z) 

-r r r . 
3( x, y, z) 


J = X| 


,-1 


J = j|Xn| 


( 18 ) 


In the _above equations superscript or subscript r refers to the reference coordi- 
nates, J is the area metric, and n is the unit surface normal. Note that all 
quantities in the above integrals are functions of the reference coordinates. Also 
for oriented bodies such as bars and beams, J and |x| may be different from each 
other if we use volume integrals throughout the sensitivity analysis. This can be 
observed in the examples discussed later in the paper. 


5. ADJOINT STRUCTURE APPROACH FOR 
GENERAL DESIGN SENSITIVITY ANALYSIS 


Discrete form 
researchers . 1 * 3 


of the adjoint variable method has been discussed by several 
Variational form of the approach based on material derivative 
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concept is described in Ref. 13 where sensitivity with respect to shape variations 
is also considered. Adjoint structure approach is described in Refs. 14-17. The 
approach has been applied to some nonlinear and shape variation problems in Refs. 

P 1 

18-20. Recently, Belegundu has traced roots of the adjoint variable method to 
methods of sensitivity analysis in optimal control literature. In addition, he has 
shown that sensitivity analysis methods for static, dynamic, shape and distributed 
parameter problems can be viewed as the general Lagrange multiplier method. This 
shows that the adjoint variable is also a Lagrange multiplier for the state 

equations which gives a sensitivity interpretation for it. This interpretation is 
extremely useful and leads to some insights into the adjoint variable method. It 
also has implications in practical applications and numerical implementations of the 
method . 

In the following derivation we combine the adjoint structure approach with the 
fixed reference volume concept to develop a general theory of design sensitivity 
analysis of linear or nonlinear structures. To avoid confusion, we use 6 and 6 to 
indicate arbitrary variations of the state fields and variations with respect to 
design variable, respectively. Also, the notation G >s will be used to indicate 

partial derivative of G with respect to ^S. Note that design sensitivity analysis 

is performed at the final state of the system denoted by left superscript t on 
various variables. Thus the virtual work equation (12) holds for the deformed 
configuration. 

Now taking variation of the functional ip in Eq. (17) with respect to design, we 
obtain 


= | 6G J dV + J G 6J dV + / 6g J df u 


+ J g S3 d? + j h 6 J dr T + / 6h J dr T + 6vp x (19) 

where <5^ represents implicit design variation of ip given as 
«1>J = /(G >s .SqS + G >e .«Q€ + G jU .6 fc u) J dV 

+ J g r «Jr d? u + I h ,„*« tu5 d *T (20) 

The basic idea of the adjoint structure approach is to replace the implicit design 
variations of the state fields in Eq. (20) by explicit design variations and the 
adjoint state fields. To accomplish this we write design variations of various 
equations as follows: 

Design Variations of the Constitutive Law (Eq. 3): 



♦,e‘*0 e + ** 


( 21 ) 
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Design Variations of Strains: 


^ E 

where 

*0 £ 


- t =t 

5 o E + 6 o £ 


( 22 ) 


|[x T ( r V6 t u T ) + ( r V«V) T X 


CT, „;t.T 


t TJ- . -T, „t T. 


+ x*( V6V) ( vV)‘x + X‘( vV) ( vSWx] 


( 23 ) 


! o e 


l[6X T ( p V t u T ) + ( r V t u T ) T 6X 


+ S x r ( r v fc u T ) ( r vV) T X + X T ( r V t u T ) ( r vV) T ax] ( 24 ) 

Here 6 represents implicit design variations of the displacements and 6 the 
explicit design variations of the strain fields. 


Design Variations of Equilibrium Equation (12): 


J SqS.c 3 J dV + / £s.Ie a J dV + / £s.e a 6J dV - J sjf.u 3 J dV 
- / Qf.u 3 6J dV - J 6qT.u 3 j dr T - | qT.u 3 6J df T = 0 


(25) 


where arbitrary variations of the primary state fields in Eq. (12) have been 
replaced by the corr esponding fields for the adjoint structure denoted by the 
superscript 'a*. The adjoint structure and the corresponding state fields are 

defined later. Substitute for 6^S from Eq. (21) into Eq. (25), use Eq. (22) and 
collect terms: 


/[ 


* 5 o e) * o s -^ a - «5 f -“ a 


+ e a .6+] J dV 


- 0 s - 6J dV - J(«o T - ua d + 5 T -“ a d ? T = 0 


Now, let us define the adjoint structure as follows: 
Loads and Boundary Conditions: 


Initial strain field 
Initial stress field 
Body force 

Specified Traction 
Specified Displacements 



G 

G 


S 

e 


f a = G 


u 


1,80 ■ \« r T 
u a0 - -g fT on r, 


(26) 


( 27 ) 
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Constitutive Law (Linear): 

S a = .(e a - e ai ) _ S ai ; S a = the adjoint stress field. (28) 

» E 

Virtual Work Equation: 

j S a .6je J dV - ( f a .6 fc u J dV - ( T a .6 fc u 3 d? T = 0 (29) 

Substitute Eq. (28) into Eq. (29): 

}(e a .* >e .6^ - e ai .* >e .6^ - S ai . 5 q£ - f a .6 fc u) J dV 

-j T a .6 fc u J d? T = 0 (30) 

Strain Field (Linear in u a ): 

e 3 = l[x T ( r Vu aT ) + ( r Vu aT ) T X + X T ( p Vu a ) (^Vfx 

+ ( p Vu a ) T X] (3D 

Substitute the adjoint equilibrium equation (30) into Eq . (26): 

* 5 S -“ a - ? 0 f -" a * cal -*,e‘ 5 0 e * 

+ S ai . + f a .6 fc u] J dV - J(^f.u a - qS. e 3 ) 6J dV 

+ J(T a .6 fc u - SqT.u 3 ) 3 d? T - J qT . u a 6J d? T =0 (32) 

Note that the variations of the state fields in Eq. (30) are arbitrary. So, 
they have been replaced as S^e = 6^e and 6^u = 6^u. 

Substitute the adjoint loads from Eq. (27) into Eq. (20): 

= J(e ai .6QS + S ai . J dV 

+ J - u a0 . 6 qT J d? u + J T a0 .6 t u3 d? T (33) 

Substitute for 5^3 from Eq. (21) into Eq. (33): 
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( 34 ) 


5^ = [(e al .<P^ £ .6 qC + e ai .<$$ + S ai .6ge + f^^ii) J dV 

+ J ~u a0 ,6 qT j dr u + / T a0 .6 t uJ d? T 

Substitute for from Eq . (22) into Eq . (3H): 

«*! = /[(e ai -* f€ + S ai ).«Q E + (e ai .* >e+ S ai ).6Q £ + e ai .«<|> 

+ f a .6 fc u] J dV + J~u a0 .5qT J d? u + j T a0 .6 T u J dr T (35) 

Substitute Eq. (32) into Eq. (35) and use Eq . (28) to obtain 
= J[«o f * ua -(e a -e al ).6+ " S a jJe - Js.Se 3 ] J dV 

+ /(gf.u 3 - qS. e 3 )<5 J dV + /(Jt.u 3 6J + sJt.u 3 J)df T 

- / s‘T .« a0 J df u (36) 


Substi tute Eq . (36) intoEq. (19): 

Sip = /[fi^.u 3 -(e 3 -e ai ).5^ - S 3 .6 qE - qS.Ie 3 + <$g]j d V 

+ /^0 f ' ua “ o S - £a + G ^ J ^ + /(*** + 8^ “ So 7 '" 30 d)d ^u 

+ J[(h + gT.u a )<$J + ( 6 h + 6 qT.u 3 ) j]d? T ( 37 ) 

Equation (37) is a general design sensitivity formula for linear and nonlinear 
structures (geometric and material nonlinearities), and shape, non-shape and 
material selection problems. Formula also gives sensitivity interpretations of 
the adjoint state fields. For example, it shows that the adjoint displacement 
field is sensitivity of functional with respect to variations of the body 
force and surface tractions. This i nterpretation has been also derived in Refs. 
21 and 22 for linear systems using the Lagrangian approach. Formula ( 37 ) also 
shows that the adjoint strain field gives variations of the functional with 
respect to the constitutive law, the adjoint stress field is related to vari- 
ations of i(i with respect to explicit design variations of the strain field, and 
variations of i|> with respect to variations of J can be recovered using adjoint 
and primary fields. These sensitivity interpretations will be observed in the 
example problems solved in the next section. These interpretations can be 
invaluable in practical applications and numerical implementations. 
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6. EXAMPLE PROBLEMS 


Several analytical linear and nonlinear examples are solved to show use of Eq. 
(37) and interpretation of various terms. Although these examples are simple, they 
can be valuable in gaining insights into numerical implementation for larger complex 
problems. Also in using Eq. (37), we will use standard symbols o for stress and e 
for strain. 


Example 1. Bar Under Self Weight 

This example is taken from Ref. 7 where sensitivity of tip displacement with 
respect to length L is calculated. We will calculate sensitivities with respect to 
all parameters of the problem to demonstrate use of formula (37) for material, 
cross-sectional and length variations. The problem definition and various trans- 
formations are shown in Fig. 3. Small displacements and linear stress-strai n law 
are assumed. The displacement field for the bar is given as u(x) = fx(2L-x)/2E 
where f is the body force per unit volume. Thus 

u (L) = fL 2 /E; 6u(L) = (L 2 /2E)6f + ( fL /E )<$L (38) 

+ ( 0) 6A - (fL 2 /2E 2 )6E 

There are at least two interpretations of this problem and both can be treated using 
Eq. (37). 

First Interpretation. In this case, Eq . (37) can be interpreted as a line 
integral with x as the only independent variable. The stress-strain law of Eq. (3) 
must be interpreted as force-strain law, as the structure is only a line element. 
Note that this must be done with the formulas given in Refs. 1A, 16, 18 and 20 when 


1 — 


t-LX. 

a = AA 

x = £L 

0 ^ 

i 

i 

t 

f 

i 

| 

-Z- 

E 

x = £L 

" J = AL 


A 

J = L 

VOLUME 
S MAPPING 

L — 

1 

1 

i 


LINE 1 

MAPPING 


1 




DESIGN VARIABLES: f, E, A, L 


Figure 3. Bar Under Self Weight 
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variations with respect to the cross-sectional area are needed. While using Eq . 
(37), the tip displacement can be treated as a boundary term or the interior term. 
We will use the latter approach. The functional for sensitivity analysis is given 
as 

* = Jl u(£)J _1 6U-1)JdC; G = uCO J _1 6(C-1 ) ; G = J _1 6(^-1) (39) 

U > U 

A 

where 5(£-1) is the Dirac delta function. The primary and adjoint fields can be 
obtained as 


uU) = fL 2 £(2-£)/2E; 

eU) = fL 2 (1-5)/E; 

e = e(5)J _1 = fL(1-C)/E; 

N = EAe = fAL(1-5); 
where N is the axial force and <|> ; 

5 ip = Jo(<5fu a - e a 5<|) - N a 


u a (^) = L£/EA 
e a (5) = L/EA 
e a = e a (5)J _1 = 1/EA 
N a = EAe a = 1 

EAe. Equation (37) reduces to 
e - N6e a + 6G)Jd5 


(40) 


+ j^fu 3 - Ne a + G)6JdC (41 ) 

Note that since we are using line integrals, the body force f = f A must be used. 
Various quantities for use in Eq . (41) are 


6<f> = ( A5E + E5A)J -1 fL 2 (1-0/E; <$f = A6f + f5A 


5G = -u(?)J 2 6(C-D6L; 


6e = e(S)6J 1 = -f(1-£)6L/E 


6e a = e a (5)5J 1 = -6L/EAL 


(42) 


Substituting all the quantities in Eq. (41) and carrying out the integrations, we 
obtain the required sensitivity equation which is the same as Eq. (38). The 
sensitivity interpretations of the adjoint fields can be directly observed. 

Second Intepretation. In this case, Eq . (3) will be treated as a volume 

integral. The functional for sensitivity analysis is given as 

4- = Jj (AL) -1 uU) 6U-1)J dAd£; G = ( AL) _1 u (5 )6 ( S"1 ) (43) 

U A 

The displacement and strain fields are the same as given in Eq . (40). However, the 
stress-strain law is the usual Hooke's Law: 
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0 = Ee = fL(1-5); 


a 3 = 1/A 


Equation (37) reduces to 


(44) 


6* = JqJ_ (6fu 3 - e a 6<j) - o a 6e - o5e a + 6G)jdAd£ 
U A 

+ jlj (fu a - oe a + G)6JdAdC 
U A 


(45) 


Various quantities for use in Eq. (45) are 

6<t> = fL6E(1-5)/E} SJ = L6A + A6L ; 6G = -(a _1 + L _1 )g (46) 

Substituting various quantities from Eqs. (40), (42) and (46) into Eq. (45), we 
again obtain the sensitivity expression given in Eq. (38). The sensitivity inter- 
pretation of the adjoint fields can be easily observed. 


Example 2. Cantilever Beam 

This example is also discussed in Ref. 7 where sensitivity of tip deflection 
with respect to the length is given. Figure 4 defines the problem and the transfor- 
mations to the reference volume. The design variables are chosen as b = (E,s,h,L). 

The tip deflection using small displacements beam theory is given as w(L) = PL^/3EI 
and its variation with respect to the design variables is given as 


6w( L) = - 


PL 


3E I 


3 3 3 2 2 

r p plV plW - pl^ - 

2 t 6E 77^2 63 77~r 6h + eT 6L 


36EI 


12EI 


(47) 




E,I,L 



V = AL = Lhs 


1 

3 — x 


x=L$ 
y = si? 
z = hC 


J = AL 


^1 




_(1 1 ) 

v 2’2' 


c 

v=i 


3- 


(1 -) 

'2 ’ 2 ) 


Figure 4. Cantilever Beam 
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The functional for the tip deflection and the function G are given as 


ip = HJ (AL) VOSU-D ALdAdC 
A 

G = (AL) -1 w(C)6U-1); G = (AL) _1 6(5-1) 

9 W 

The primary and adjoint structure solutions are given 

W ,££ = EI~ ^ 1- ^ ; w ^) = 6 El d'Z') 

3 r 2 


W ,£E; El w 6EI (3 

The sensitivity formula of Eq . (37) is reduced to 

= Jq}_ (-e a 64> - ale 3 - o 3 6e + 6G)ALdAd£ 


+ Inl_ ("oe 3 + G)6(AL) dAd^ 

u A 


(48) 

(49) 


(50) 

(51) 


(52) 


The following quantities are needed to complete integrations in Eq . (52): 
e = £w hL 2 ; e a = £w 3 hL 2 ; SG = w5 ( £- 1 ) <5 ( AL) 1 

» S>S> > SS) 

a = Ee; 6(J> = eSE; o 3 = Ee 3 

6e = £w 6(hL _2 ); 6e 3 = £w 3 J(hL~ 2 ) 

» 9 SS 

Substituting these quantities in Eq. (52) and carrying out the integrations we get 
the sensitivity expression given in Eq . (47). It is interesting to again note that 
the adjoint displacement field given in Eq. (51) represents the sensitivity of the 
primary displacement field (Eq. 52) with respect to the load parameters P; i.e. 

u a (5) = dip/ dP. 


Example 3. Materially Nonlinear Problem 

Consider the bar of Fig. 3 subjected to a load P in the x direction at the free 

end. The material for the bar obeys a nonlinear stress-strain law o = Ee^ 2 (e>0); 

1 /2 

so (j> = Ee . We will consider E, A and L as design variables and determine sensi- 
tivity of the tip deflection. Transformation to the reference volume gives x = L£, 
a = AA, J = AL, V = AL, and V = 1. Nonlinear analysis of the primary structure 
y i el ds : 
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»«> ■ 4l- ' u(L > - - 2 4k * ■ z 4i «* * 4~, «L 
A E A E J A^E A E 

The functional for sensitivity analysis is given as 

* - JJJjA" 1 l” 1 u( 5)«(5-1) (AL)dAdc 

G = a _1 l _1 u(5)6(c-i); G = A _1 l" 1 8(5-1) 

» n 

The adjoint structure is linear with the stress-strain law as 


a 1 _ 1/2 a 
i e = r Ee e 
,e 2 


AE 2 a 
2P E 


The equilibrium equation for the adjoint structure gives 

a 2PL 


a,_\ 2PL£ 
u ( 5 ) = 

AT; 


,5 


2 2 
A E 


(53) 

(54) 

(55) 

(56) 

(57) 


The sensitivity formula of Eq. (37) reduces to 

8 * = JqJ_ C~ e a 64 > ~ o 3 le - ale 3 + 6 G)jdAdC 


Jq|_ (~ ae 3 + G)sJdad£ 


(58) 


Various quantities for Eq. (58) are 


e = u L 

9 5 


-1 PL 


2 2 ’ 
AE 


e 3 = u 3 L 1 = 6 G = u8U-1)6(AL) 1 

’ * AE 


o = Ee 1/2 ; <p = (1/2)Ee 1/2 , S* = e 1/2 6E; a 3 = <j> e 3 
f € , e 


6 e = u 6 L _1 = -L 2 u 6 L; le 3 = -L 2 u 3 6 L 

9 S 9 S 9 

Substituting these quantities into Eq . (58), we obtain the sensitivity formula of 
Eq . (53). It is interesting to observe sensitivity meaning of the adjoint displace- 
ment field in Eq. (57); i.e. u 3 (£) = d\p/ dP . 


Example 4. Geometrically Nonlinear Problem 

Consider the two bar structure shown in Fig. 5. The material for the structure 
is linear, so a = Ee. Transformation to the reference volume is shown in the 


1 27 


figure. The design variables for the problem are b = ( E, A, L) . The strain for the 
problem is given as 

t ft. o. i . 0 , ‘1 2 , -2 

Q e = [ L - LJ/ L = - w L » e ( 59 ) 

The deflection at the center and member strains are calculated as 
P 1/3 L - 2 12 P 2 / 3 L 2 

W jjz; e - EL = - W = yjr (60) 

(EA) ' /3 2 (EA ) 2/3 

The incremental equilibrium equation in terms of displacement at the center is 
2 -3 

3EAw L 6 w = 5P. The functional for sensitivity analysis is given as 


* = ( AL ) 1 w(5)6(5-1)ALdAdC 


G = (AL) Vo-SU-D; G , = (AL) 1 5 ( C _ 1 ) 

« W 


The equilibrium equation for the adjoint structure (using the incremental equili- 
brium equation of the primary structure) is given as 


3 EAL Vw a = |qJ_(AL) 1 6 ( 1 ) ALdAd£ = 


3EAw 2 3E 1 /3 A 1 /3 P 2/3 

Total axial displacement and displacement at any point are given as 


u a (L) = w a sine = w a wL 1 ; u a (£) = w 3 wL \ 


P 



a = A A 


Figure 5. Two Member Structure 


The adjoint strains are given as 

a a . -2 -a. -2 -a a 

e = w wL = el , e = w w 

The sensitivity formula of Eq. (37) reduces to 

= 2 (-e a 64> - o a 5e - ale 3 + 6G)jdAdC 
0 A 

+ 2jH_(-ae 3 + G)6JdAd5 

^ A 


( 65 ) 


( 66 ) 


Note that factor of 2 is used because volume integrals in Eq. (37) are for the 
entire structure. Various quantities for Eq . (66) are 

5<J> = e6E ; le = e6L~ 2 = -w 2 l“ 3 6L; le 3 = e 3 6L~ 2 = -2 ww 3 L _3 6L 
0 = Ee = ^Ew 2 L -2 ; a 3 = Ee 3 = Eww^" 2 ; 6J = L6A + A6L (67) 


Substituting these quantities into Eq. (66), we get 


6i|> 


r u h 

3 E 1/3 a 1/3 


«E - 


p 1/3 l 

3E 1/3 A ,|/3 


5A + 


(P/EA) 1/3 Il 


( 68 ) 


which can also be obtained directly from Eq . (60). Comparing w and w a in Eqs. (60) 
and (63), we again observe the sensitivity interpretation of the adjoint displace- 
ment field. 

2/3 1 /3 2/3 

Next, consider the member stress given in Eq . (67) as a = (P E )/(2A ). 

Its design sensitivity is given as 


p 2/3 _ p 2/3 E 1/3 _ 

6a = -jjt-J To « e - ' -Tk- 6A+(0)SL 

6E 2/3 A 2/3 3A 5/3 


(69) 


The functional for design sensitivity analysis is given as 


4) = JqJ_(AL) _ 1 a AL dAdS = /q}_G(<T)AL dAd^ 


(70) 


The adjoint load G, u in this case is zero but initial strain in the adjoint 

3.1 — 1 q 331 

structure and stress strain law are given as e = G = (AL) , o = E(e -e ). 

a 

The adjoint equilibrium equation in terms of central displacement gives w 

(wA) 1 L / 3 and e = ww L = (AL) /3. Substituting appropriate quantities in Eq. 
(37), it can be verified that Eq . (69) is obtained. It can be also directly 

verified that w 3 = dip/ dP = da/dP. 


129 


7. CONCLUDING REMARKS 


A general formula for design sensitivity analysis of linear and nonlinear 
structures using variational approach has been developed. Equations of continuum 
mechanics are used and the concepts of reference volume and adjoint structure are 
exploited. Use of the formula is demonstrated on a few simple analytical problems. 
The theory can be easily adapted for finite element modelling of structures. The 
finite element models for the primary and adjoint structures can be independent of 
each other. For modelling of design optimization problems, the concept of a refer- 
ence volume is translated to the concept of a design element that is invariant with 
respect to design changes. These observations can have considerable implications in 
numerical implementations for design sensitivity analysis and optimization of 
complex structures. 

Considerable numerical work has been done for design sensitivity analysis and 

optimization of linear structures.^ Material derivative approach has been exploited 

for shape optimization. In this regard recent work of Choi and Co-workers , 2 ^ Yang 

and Co-workers 24 and Hou and Co-workers 25 is significant. Yang and Co-workers 24 
have shown equivalence of variational and finite element formulations of design 
sensitivity analysis of shape problems for linear structures. This equivalence can 

also be shown for nonlinear problems. Hou and Co-workers 25 have discussed some 
difficulties with the material derivative approach of design sensitivity analysis of 
linear shape problems. They have suggested numerical procedures to improve accuracy 
of the approach. 

Design sensitivity analysis and optimization with nonlinear response is just 
beginning to be studied. Finite element approach for nonlinear stresses, strains, 

displacements and the buckling load has been recently studied. 2 ^ ^ 2 . More research 
needs to be done to fully develop this area. 
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